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We investigate the Unruh effect on entanglement taking into account the spin degree of freedom of 
the Dirac field. We analyze spin Bell states in this setting, obtaining their entanglement dependence 
on the acceleration of one of the partners. Then, we consider simple analogs to the occupation num- 
ber entangled state |00) + |11) but with spin quantum numbers for |11). We show that, despite their 
apparent similitude, while the spinless case is always qubitxqubit, for the spin case acceleration 
produces a qubitxqu4it state. We also introduce a procedure to consistently erase the spin informa- 
tion from our setting preserving occupation numbers. We show how the maximally entangled state 
for occupation number emerges from our setting. We as well analyze its entanglement dependence 
on acceleration, obtaining a greater entanglement degradation than in the spinless case. 
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I. INTRODUCTION 

Despite their apparently separated application areas, 
general relativity and quantum information are not dis- 
joint research fields. On the contrary, following the pio- 
neering work of Alsing and Milburn [1] a wealth of works 
0-13, fflHI has considered different situations in which 
entanglement was studied in a general relativistic setting, 
for instance, quantum information tasks in the proximity 
of black holes_M H, [m [13 ! entanglement in an expanding 
universe [71, |l3| , entanglement with non-inertial partners 
I, Hi etc. 

Entanglement behavior in non-inertial frames was first 
considered in [l| where the fidelity of teleportation be- 
tween relative accelerated partners was analyzed. Af- 
ter this, occupation number entanglement degradation of 
scalar Q and Dirac [6i] fields due to Unruh effect J, 4, [l^ 
was shown. Recent works studied the effect of the instan- 
taneous Wigner rotations and Thomas spin precession on 
entanglement [ll] , 17 1 . 

The previous work on Unruh effect for Dirac field 
mode entanglement does not consider the spin of the par- 
ties. Hence, only two occupation numbers n = (0, 1) are 
allowed for each mode. Higher values of n are forbidden 
by Pauli exclusion principle. However, addressing the ef- 
fect of Unruh decoherence on spin entanglement can only 
be done by incorporating the spin of the parties in the 
framework from the very beginning. As a consequence, 
occupation number n = 2 is also allowed. This fact will 
affect occupation number entanglement which has to be 
reconsidered in this new setting. For this purpose, in this 
work we will study the case of two parties (Alice and Rob) 
sharing a general superposition of Dirac vacuum and all 
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the possible one particle spin states for both Alice and 
Rob. Alice is in an inertial frame while Rob undergoes a 
constant acceleration a. 

We will show that Rob -when he is accelerated re- 
spect to an inertial observer of the Dirac vacuum- would 
observe a thermal distribution of fermionic spin 1/2 par- 
ticles due to Unruh effect [13] ■ Next, we will consider 
that Alice and Rob share spin Bell states in a Minkowski 
frame. Then, the case in which Alice and Rob share 
a superposition of the Dirac vacuum and a specific one 
particle state in a maximally entangled combination. In 
both cases we analyze the entanglement and mutual in- 
formation in terms of Rob's acceleration a. 

Finally, we will study the case when the information 
about spin is erased from our setting by partial trac- 
ing, remaining only the occupation number information. 
Here, entanglement is more degraded than in This 
comes about because more accessible levels of occupation 
number are allowed, so our system has a broader margin 
to become degraded. 

This paper is structured in the following sections. 

In sectionsini and lllll we introduce the basic formalism 
and notation to deal with Dirac fields from the point of 
view of an accelerated observer taking its spin structure 
into account. In section ITVl we study how the Minkowski 
vacuum state is expressed by an accelerated observer 
when the spin of each mode is included in the setting, 
discussing the implications of the single-mode approxi- 
mation often carried out in the literature Also, we 
build the one particle state with spin s in Rindler coordi- 
nates and analyze the Unruh effect when the spin struc- 
ture is included. Here we discuss the necessity of trac- 
ing over Rindler's region IV for Rob, since it is causally 
disconnected from region I in which we consider Rob's 
location. In section |V] we analyze how entanglement of 
spin Bell states is degraded due to Unruh effect. We show 
that, even in the limit of a — >■ oo, some degree of entangle- 
ment is preserved due to Pauli exclusion principle. Then 
we analyze Unruh effect on a completely different class of 
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maximally entangled states (like |00) + \ss') where s and 
s' are z component of spin labels) comparing it with the 
spin Bell states. In section IVTl we show that the erasure 
of spin information, in order to investigate occupation 
number entanglement alone, requires considering total 
spin states for the bipartite system. Finally, our results 
and conclusions are summarized in section I VIII 

II. THE SETTING 

We consider a free Dirac field in a Minkowski frame 
expanded in terms of the positive (particle) and the neg- 
ative (antiparticle) energy solutions of Dirac equation no- 
tated V'/i^^ and ^p^ ^ respectively: 

= ^ /'d3fc(a,,,^+^+6^^-J (1) 

s •' 

Here, the subscript k notates momentum which labels the 
modes of the same energy and s = {t, i} is the spin label 
that indicates spin-up or spin-down along the quantiza- 
tion axis, ak^s and bk,s are respectively the annihilation 
operators for particles and antiparticles, and satisfy the 
usual anticommutation relations. 

For each mode of frequency k and spin s the positive 
and negative energy modes have the form 

V'^, = -i=u^(fe)e±^('=— '^o*) (2) 

where uf{k) is a spinor satisfying the normalization re- 
lations ±uf{k)uf,(k) = {ko / 171)6 ss',uf{k)uf,{k) — 0. 

The modes are classified as particle or antiparticle re- 
spect to dt (Minkowski Killing vector directed to the fu- 
ture). The Minkowski vacuum state is defined by the 
tensor product of each frequency mode vacuum 

\0)=^\0k)^\0k')- (3) 

k,k' 

such that it is annihilated by ak.s and bk.s for all values 
of s. 

We will use the same notation as reference where the 
mode label will be a subscript inside the ket, and the ab- 
sence of subscript outside the ket indicates a Minkowski 
Fock state. 

In this way, and as a difference with previous works, 
we will consider the spin structure for each mode, and 
hence, the maximum occupation number is two. This 
introduces the following notation 

4A,s'\0) = \s4)Ss.-s' (4) 

If s = s' the two particles state is not allowed due to Pauli 
exclusion principle, so our allowed Minkowski states for 
each mode of particle/ antiparticle are 

{\Okf ,\tkf Aikf AUkf}. (5) 



Consider that we have the following Minkowski bipar- 
tite state 

\A.,k,) = /i|0fej+|0fcj+ + a|t.j+|tfcj++/3|tfej + 

(6) 

with fi = v^l - |a|2 - |/3|2 - |7|2 _ The subscripts 
A, R indicate the modes associated with Alice and Rob 
respectively. All other modes of the field are unoccupied 
-that is to say that the complete state would be |<I>) = 

\4'k^,ka) ^ |Ofe)^ |Ofc)"]-- 

This state generalizes the Bell spin states (for instance, 
we have | (/)"'") choosing a = S = l/'v/S) or a modes entan- 
gled state (for instance choosing a = ^ = 1/V2). With 
this state ([6]) we will be able to deal with two different 
and interesting problems at once, 1. Studying the Unruh 
decoherence of spin entangled states and 2. Investigat- 
ing the impact of considering the spin structure of the 
fermion on the occupation number entanglement and its 
Unruh decoherence. 

Later on, under the single mode approximation, we 
will assume that Alice is stationary and has a detector 
sensitive only to the mode kA, and Rob moves with uni- 
form acceleration a taking with him a detector sensitive 
to the mode kR. 

III. RINDLER METRIC AND BOGOLIUBOV 
COEFFICIENTS FOR DIRAC FIELDS 

An uniformly accelerated observer viewpoint is de- 
scribed by means of the well-known Rindler coordinates 
p^ . In order to cover the whole Minkowski space-time, 
two different set of coordinates are necessary. These sets 
of coordinates define two causally disconnected regions 
in Rindler space-time. If we consider that the uniform 
acceleration a lies on the z axis, the new Rindler coordi- 
nates {t, X, y, z) as a function of Minkowski coordinates 
{i,x,y,z) are 

at = e"^ sinh(ai), az = e"^ cosh(ai), x — x, y = y (7) 

for region I, and 

at ~ — e"^ sinh(at), az — — cosh(at), x = x, y = y 

(8) 

for region IV. As we can see from fig. 1, although we have 
covered the whole Minkowski space-time with these sets 
of coordinates, there are two more regions labeled II and 
III. To map them we would need to switch cosh o sinh 
in equations (O,®. In these regions i is a spacelike co- 
ordinate and z is a timelike coordinate. However, the so- 
lutions of Dirac equation in such regions are not required 
to discuss entanglement between Alice and an accelerated 
observer since he would be constrained to either region I 
or IV, having no possible access to the oppos ite regions 
as they are causally disconnected [1, 0, [l9l. |20|| . 
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IV 




III 

Figure 1: Rindler space-time diagram: lines of constant posi- 
tion z — const, are hyperbolae and all the curves of constant 
proper time t for the accelerated observer are straight lines 
that come from the origin. An uniformly accelerated observer 
Rob travels along a hyperbola constrained to region I 

The Rindler coordinates z,t go from — oo to oo in- 
dependently in regions I and IV. It means that each 
region admits a separate quantization procedure with 
their corresponding positive and negative energy solu- 
tionsi and {C^,^"}- 

Particles and antiparticles will be classified with re- 
spect to the future-directed timelike Killing vector in each 
region. In region I the future-directed Killing vector is 

d't = + l^a^ = a{Sdi + id,), (9) 
whereas in region IV the future-directed Killing vector is 

or = -di. 

This means that solutions in region I, having time de- 
pendence i/)^"'' ^ g-^kot -> 0^ represent positive 
energy solutions, whereas solutions in region IV, having 
time dependence ipl'^ ~ g-»fcot ^j^j^ -> ^re actually 
negative energy solutions since d^^ points to the opposite 
direction of df [1, HO] • As I and IV are causally discon- 
nected ipk^s^ ^^'^ ^fe^ ^'^^y have support in their own 
regions, vanishing outside them. 

Let us denote {ci,k,s,c\ ). ^) the particle annihilation 

and creation operators in region I and {di,k,s,dj 
the corresponding antiparticle operators. Analo- 
gously we define {civ,k,s,cly,. ,,,div,k,s,d\y f. ,,) the par- 
ticle/antiparticle operators in region IV. 

These operators satisfy the usual anticommutation re- 
lations {cR,k,s,cli',k',s'} = SjiR'Skk'Sss' where the sub- 



^ Throughout this work we will consider that the spin of each mode 
is in the acceleration direction and, hence, spin will not undergo 
Thomas precession due to instant Wigner rotations [2ll| . 



script R notates the Rindler region of the operator R = 
{I, IV}. All other anticommutators are zero. That in- 
cludes the anticommutators between operators in differ- 
ent regions of the Rindler space-time. 

Taking this into account we can expand the Dirac field 
in Rindler coordinates analogously to (IT|): 

^ ^ J2 [d'k {ci,k,s^i++4 ,^ ji-+civ.k.si'Z'' 

+d\v.k.M^-) ■ (10) 

Equations (P) and ((TU)l represent the decomposition of 
the Dirac field in its modes in Minkowski and Rindler 
coordinates respectively. We can relate Minkowski and 
Rindler creation and annihilation operators by taking ap- 
propriate inner products [1, [20l - |22j . The relationship be- 
tween Minkowski and Rindler particle / antiparticle oper- 
ators is linear and the coefficients that relate them are 
called Bogoliubov coefficients: 

ak,s = cosrci^k,s-e'''^smrd\y_^._^ 

bL = cosr4^,,, + e-*'^sinrc/,-fe,-, (11) 

where 

tanr = e"''-^ (12) 

and is a phase factor that will turn out to be irrelevant 
for our purposes. Notice that as we are working with 
two spatial-temporal dimensions and with massless Dirac 
field, the relation between Rindler modes and Minkowski 
modes is given in ([TT|) . We will discuss in the conclu- 
sions the implications of considering extra dimensions 
and massive fields, where Minkowski modes are spread 
over all positive Rindler frequencies . 

Notice from Bogoliubov transformations dTTI) that the 
Minkowski particle annihilator ak,s transforms into a 
Rindler particle annihilator of momentum k and spin s 
in region I and an antiparticle creator of momentum —k 
and spin — s in region IV (in region IV all the magnitudes 
that arc not invariant under time reversal change). 

IV. UNRUH EFFECT FOR FERMION FIELDS 
OF SPIN 1/2 

Now that we have the relationships between the 
creation and annihilation operators in Minkowski and 
Rindler coordinates, we can obtain the expression of the 
Minkowski vacuum state for each mode |Ofe) in Rindler 
coordinates. For notation simplicity, we will drop the k 
label in operators/states when it does not give any rel- 
evant information, but we will continue writing the spin 
label. 

It is useful to introduce some notation for our states. 
We will denote with a subscript outside the kets if the 
mode state is referred to region I or IV of the Rindler 
space-time. The absence of subscript outside the ket 
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will denote Minkowski coordinates. The ± label of 
particle/antiparticle will be omitted throughout the pa- 
per because, for the cases considered, a ket referred to 
Minkowski space-time or Rindler's region I will always 
denote particle states and a ket referred to region IV will 
always notate antiparticle states. 

Inside the ket we will write the spin state of the modes 
as follows 



\S)I = CI\0}J, \s)jy=d\yJO} 



IV 



(13) 



which will notate a particle state in region I and an an- 
tiparticle state in region IV respectively, both with spin 
s. 

We will use the following definitions for our kcts 

\n)j - 4Aio), = -4,4^10), 

\n)lV = 4y^d\y^\0),y=-d\y/,y^\0},y (l4) 

and, being consistent with the different Rindler regions 
operators anticommutation relations, 

k)/ W)iv = 44y.' |0)/ 10),,, = -4ys'4 |0)/ |0),^ 



4vs' k)/ |0)/y = - W)jy 



(15) 



Now, it is useful to note that (|lll) could be expressed 
as two-modes squeezing transformation for each fc [J, Q 



as.k 



k.s 



s 



where 



exp 



{'I 



''IV, -k, 



Cl.k.s 



sC -I- CTk^ 



■ "'IV,- 



(16) 



J't> 



(17) 

So, analogously to [1, @|, it is reasonable to postulate 
that the Minkowski vacuum is a Rindler two-mode parti- 
cles/ antiparticles squeezed state with opposite spin and 
momentum states in I and IV. Contrarily to Iq, con- 
sidering that the modes have spin, occupation number 
for each k is allowed to be 2, being higher occupation 
numbers forbidden by Pauli exclusion principle. In the 
literature the analysis is restricted only to one mode of 
the Minkowski field, but we can restrict our analysis to 
some sector of the Minkowski vacuum ([3]), defining for 
the particles sector 



10, 



(18) 



fei, 



such that the particle sector of can be written as 

|o) = |6)^(8)p^fc^,...,feJ0p). 

In this fashion we are considering a discrete number 
n of different modes fci,...,fc„, so Minkowski vacuum 
should be expressed as a squeezed state in Rindler coor- 
dinates which is an arbitrary superposition of spins and 



momenta. This will be useful to discuss what would hap- 
pen if we relax the single-mode approximation carried out 
often in the literature and let our detectors have a small 
mode spread. 

|o) = c°|o),|o),^ + ^c,\,,Ji)Ji),^ 

Sl 

+ 51 ^Si,S2,ki,k2^sl:S2 |2)jy + ■ • • 

Sl,S2 

ki ,k2 



Si,...,Sn 

k-L,...,kn 



^ / y '-"Si,...,S2n,kl,...,k2n^Si,. 



,S2n 



Sl,...,S2„ 
ki,...,k2r, 



Where, here, the notation is 



2n 



2n) 
I IV 



(19) 



|7Tl)j — l^l ;Sm jfcm) / I ^i, fci ; . . . ; S^, ^nijjy 



with 



(20) 



being 



J2 {m\iv{rh'\j\rh)j\rh)jy^im\f (22) 



Si,...,Srr 

k-L,...,k„ 



and the symbol ^ is 

^ki,...,kTn - 



If Sj = Sj and ki = kj i ^ j 

1 Otherwise 



(23) 



which imposes Pauli exclusion principle constraints on 
the state (quantum numbers of fermions cannot coin- 
cide). 

Notice a pair of aspects of this notation for the mul- 
timode case. First, in the series in all the possible 
orders of the operators are implicitly written. Due to 
the anticommutation relations of the fermionic operators, 
terms with different orderings of the creation operators 
are related, i.e. 



\si,ki;s2, fe)/ |-si, -fci; -S2, ~k2)iv 

|S2, fe; Si.kijj |-S2, -k2] -Sl, -kij^y 



(24) 



So, without loss of generality, we could choose not to 
write all the possible orderings of the operators in p^ . 
The difference between taking all the possible orderings 
of the operators into account and taking only one repre- 
sentant is a factor m! in the constants C™. From ([24]) we 
can also see that the coefficients C™ are symmetric with 
respect to s^, ki index permutations. 

Second, as there are only n different modes 
(/ci, . . . , fc„), the last summation in equation (I19p has 
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only {2n)\ terms due to ([23)) . These terms are all the 
different permutations of the creation operators for pairs 
of opposed spins for each mode. There would be only 
one summand-instead of (2n)!-in the simplified notation 
where we do not write all the different permutations of 
the operators but only one representant. It means that, 
in this simplified notation, the series of terms with n pairs 
has the same summands as the series with the vacuum 
state (i.e. only one). Actually, in this notation-i.e. if we 
count all the different order permutations as only one- 
the series with C^" to C""*"^ has exactly the same number 
of summands as the series with C° to C"~^. 

To obtain restrictions on the values of the coefficients 
C™ we demand that the Minkowski vacuum has to be 
annihilated by the particle annihilator, ak^^so |0) = 0- 
Translating this into Rindler coordinates we have 



e*"^ sin r d 



t 

IV -ko 



|0) = (25) 



where the vacuum state should be expressed in Rindler 
coordinates using 

As the elements form an orthogonal set, from ([^5)) 
we see that all the terms proportional to different ele- 
ments of the set should be zero simultaneously, which 
gives the following conditions on the coefficients 



• J, as a function of C° 



Clkocosr-C°e''''smr = 



(26) 
(27) 



since equations (I^B1) . (P7|) should be satisfied V/cqj we ob- 
tain that k — k — const, since does not depend 
on k or s. We will denote = C^. 



• Cc c i. i. as a function of 

C'e'*8mr~2Cl,^ko.kCOsr = 



since equations 



(28) 
(29) 



tain that Cl , u u 

Si ,S2 ,/C2 



should be satisfied Vfco: we ob- 
where does not depend on 
spins or momenta since does not depend on k or s, the 
only dependence of the coefficients ([T5|) with ki and st is 
given by the Pauli exclusion principle, this dependence 
comes through function (1^51) . 

In fact it is very easy to show that all the coefficients 
are independent of si and ki -apart from the Pauli ex- 
clusion principle constraint.- Using the fact that C° does 
not depend on Si and ki and noticing that by applying the 
annihilator on the vacuum state and equalling it to zero, 
we will always obtain the linear relationship between C" 
and C"~^ given below. 

• C™ as a function of C"~i 

C""-^e*'^sinr-mC""cosr = (30) 
™-ie*'^sinr-mC""cosr = (31) 



C 



We finally obtain that C" is a constant which can be 
expressed as a function of C° as 



C" = — e™'^tan™r 
m! 



(32) 



And then, vacuum state can be expressed as 

|0) = C^\0)j\Q),y + C^Y.\^)^\l)^ 



'IV 



Si 

ki 



Sl,S2 
fcl,fc2 

Sl....^Sn 
kl . . ..,krL 



IV 



n)i \n)jy + 



n), n 



IV 



(33) 



Sl,...,S2„ 

ki,....k-2„ 



where the only parameter not fixed yet is . 

We can now fix by imposing the normalization of 
the Minkowski vacuum in Rindler coordinates (0 |0) = 1, 
this condition can be written as 



0|2 



Where 



2n 



ET™tan2™r+ E T2„-™ tan^™ r 



m— Ti+1 



T - ffei,...,fe„ 
™ ~ / , Ssi,...,s„ 



= 1 

(34) 
(35) 



fci , kj. 



and we have defined To = 1. This expression gives (for- 
mally) the value of (except for a global phase) when 
considering the populated levels in an arbitrary number 
of modes of the field 



T 



tan^" r 



2n 

E 

m— n+l 



Tsn-^tan^^r 



-1/2 



(36) 

We can see that if we take the limit a— >0=>r— >Owe 
recover the Minkowski Vacuum. We will come back to 
the behavior at the limits below. 

Notice that the state ([33]) is only normalizablc in this 
discrete limit. This comes about because the Minkowski 
and Rindler representations are not unitary equiva- 
lent (there is not unitary operator connecting the two 
vacua) ^. It prevents us from taking the continuous limit 
in the expression but does not invalidate the treat- 
ment as equation (|33[) can be considered as a superpo- 
sition of a finite number of individual modes which are 
perfectly well defined 0. 



this is a old-known problem (see reference '2^ chapter 2) 
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To address the problem of showing how the presence 
of spin degrees of freedom affects the entanglement be- 
tween accelerated observers, it is useful to use the single 
mode approximation [T, 18] in the same way that in 
This is valid if we consider Rob's detector so sensitive to 
a single particle mode in region I that we can approxi- 
mate the frequency uja = ^Oa observed by Alice to be 
the same frequency observed by Rob, uia ~ ^r- As a 
consequence, the populated levels we are looking at are 
in this single frequency @ (See also discussion in [l^) 
so we can consider the sums over ki in just like a 
sum of only one mode k = uja- This is equivalent to re- 
stricting the analysis to the sector |Ofc) of the complete 
vacuum ([3]). Since the goal of this work is to show the 
effect of spin degrees of freedom on the entanglement 
for non-inertial observers, this approximation allows us 
to compare our results withprevious literature on scalar 
and spinless fermion fields [j, Q . 

We have to notice that since the observer Rob is ac- 
celerated, his possible measurements are affected by a 
Doppler-like effect. Given that the velocity of the ob- 
server is V = at = t/ z = tanh(ai), the Doppler effect will 
shift the sensitivity peak of the detector. Namely, if at 
the instant t = Q his detector is sharply tuned to a fre- 
quency uj£) — UJA, to compensate this Doppler effect at 
some instant t — t the detector should be tuned to the 
frequency ljd = c^'^wa- This implies that any detector 
will eventually become insensitive to the populated lev- 
els of the Minkowski field. In order to do the theoretical 
analysis in this work, we can consider either that his de- 
tector can be sharply tuned to the frequency e°''^LUA for 
each instant, or that Rob has a set of individual detectors 
each one sharply tuned to the proper frequency for each 
instant. 

Carrying out this single-mode approximation, the 
Minkowski vacuum for a single mode is a Rindler two- 
mode particles/antiparticles squeezed state with oppo- 
site spin states in I and IV. Considering that the modes 
have spin, occupation number is allowed to be 2 for each 
fc, being higher occupation numbers forbidden by Pauli 
exclusion principle. As a consequence 

|0) = V\0)j\0)jy+A\t)j\i)jy+B\i),n)jy 

+c\n)j\n)jv (37) 

Notice that V is the analogous to C°, A and B are anal- 
ogous to and C| respectively and C is analogous to 
in the expression but considering only one rep- 
resentant of all the 2 possible orders for the pair . 

To obtain the values of the coefficients V, A, B, C we 
demand that the Minkowski vacuum has to be annihi- 
lated by the particle annihilator, |0) = 0. Translating 
this into Rindler coordinates we have 



which implies 



= 



cos r c/.s — e*"^ sin r d\ 



IV.- 



[^|0),|0),^ 



cosr[A\Q)j\i)jy + 

+ C{5s^\\)j\n)iv 
-e"!' sinr[F|0)^|-, 

+i?<5,t i;>/ \n)iv] • 



S<5.;|0),|t),v 
-5si |t>/in>/y) 



I IV 



A5 si it>,in> 



IV 



(39) 



This equation gives 4 conditions (two for each value of 
s), although only 3 of them are independent 



A cos r — Ve^'^ sin r = 

C cos r — Be"^ sin r = 

B cos r — Ve^'^ sin r = 

Ccoar - Ae">'smr = 



A = B^ Fe*"*tanr 
C = Ve^'^ tan^ r 



(40) 

To fix V we impose the normalization relation for each 
field mode (0|0) = 1 => \V\^ = 1 - - \B\^ - |Cp, 
imposing this we finally obtain the values of the vacuum 
coefficients. 



V 
A 
B 
C 



g40 gjjj J. (,Qg ^ 



(41) 



A It)/ \i)iv + B |t)/y + C \n)i in)/y](38) 



Notice that comparing this result with expressions ([32)) 
and (1551) . as we have truncated the series in the 
value of V will be different from the case when more 
than one mode is considered -C" instead of V-. If we 
restrict the series on m to only one mode n = 1 in p4p . 
we obtain that C° — > 1/(1 + tan^ r) — cos^ r and we get 
then proper values for A = B = and C = 2!C^. 

Since from (jl2p a^oo=>r— >7r/4, comparing V 
with (j36p . we can see that, while under the single-mode 
approximation the limit of infinite acceleration leads to a 
finite distribution of the Minkowski vacuum over Rindler 
states, when considering the multimode Rindler expres- 
sion for the vacuum state (j33p the combined limit of 
n oo and a— )-oo=^r— >7r/4 leads to a complete 
fading away of the amplitudes over all the Rindler modes 
as Co — > 0, which may not be the case for finite a. This 
is beyond the scope of this article but we will discuss in 
the conclusion that it may have very strong implications 
on the entanglement of fermionic fields for accelerated 
observers. 

So finally, under the single mode approximation, the 
Minkowski vacuum state in Rindler coordinates is as fol- 
lows 

|0) = cos^r |0),|0),^+e'^sinr cosr(|t),|;),^, 

+ i;)/lt)zy) + e2'^sin2r \n) AU) iv (42) 

Now we have to build the one particle (of spin s) state 
in Rindler coordinates. It can be readily done by ap- 
plying the Minkowski particle creation operator to the 
vacuum state \s) = a| |0), and translating it into Rindler 
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coordinates: 

\s) = 



sin r djY_- 



''r\0)j\0)jy 



W'sinrcosri\t)j\i)jy + \i)jn)jv) 

+e^^'f' sin' r\n) I \n)jv] (43) 

That means 

It) = cosr|t),|0),^, + e'*smr|n),|t)/y 

i;) = cosr|;),|0),^,-e^^sinr|n>,|;),v (44) 

The three Minkowski states |0) , |t) i li) correspond to 
the particle field of mode k observed by Alice. However, 
since Rob is experiencing a uniform acceleration he will 
not be able to access to field modes in the causally dis- 
connected region IV, hence, Rob must trace over that 
inaccessible region as it is unobservable. 

Specifically, when Rob is in region I of Rindler space- 
time and Alice observes the vacuum state, Rob could 
only observe a non-pure partial state given by pa = 
Tr/y (|0) (0|) that is 

PR = cos^ r \0)j{0\ + sin^ r cos^ r (|t),(t| 

+ i;),ai) + sinVin),(ni (45) 

But while Alice would observe the vacuum state of mode 
fc, Rob would observe certain statistical distribution of 
particles. The expected value of Rob's number operator 
on the Minkowski vacuum state is given by 

(Nr) = (0 \Nr\ 0) = Tvijv {Nr |0)(0|) = Tii (Nrpr) 



= Trr 



(4^C7t + c|^c/;) PR 

Substituting the expression pS)) we obtain 

(Nr) -2sin2r 
using ([T^ we obtain that 



where A;_b is the Boltzmann's constant and 



27rfcf 



(46) 



(47) 



(48) 



(49) 



is the Unruh temperature. 

Equation (j48]) is known as the Unruh effect [13, [l!] , 
which shows that, for a two-dimensional space-time, 
an uniformly accelerated observer in region I detects a 
thermal Fermi-Dirac distribution when he observes the 
Minkowski vacuum. We obtain a factor 2 contrarily to 
Ref. [6] due to the degeneracy factor 25* + 1. 



V. SPIN ENTANGLEMENT WITH AN 
ACCELERATED PARTNER 

In previous works 0, IB| it was studied how Unruh de- 
coherence affects occupation number entanglement in bi- 
partite states as 



|M/) = -i=(|00)-f 111)) 



(50) 



where the figures inside the kets represent occupation 
number of Alice and Rob modes respectively, barring any 
reference to the spin of the field modes. 

Here, where we have included the spin structure of each 
mode in our setting from the very beginning, it is possible 
to study the effects of acceleration in spin entanglement 
decoherence, which is different from the mere occupation 
number entanglement. 

First of all, we build a general bipartite state that could 
be somehow analogous to state (|50p studied in [6] , limit- 
ing the occupation number to 1 but including the spins 
of each mode. 

I*) = M|OA)|Ofl.)+a|tA)|t7?)+/3|tA)|;ii) 

+7 \U) Itfi) + 5 \U) \iR) (51) 



with p = — lap 



171' 



I (5 p. The subscripts 



A, R indicate the modes associated with Alice and Rob 
respectively. We will suppress the labels A, R from now 
on, and we will understand that the first character in a 
ket or a bra corresponds to Alice and the second to Rob: 
\s,s')^\sa)\s'r) 

This general setting ([5T|) allows us to study in this sec- 
tion what happens with spin entanglement under acceler- 
ation of Rob and also what happens with the occupation 
number entanglement when considering sates analogous 
to ((50)) but taking the spin structure into account. It will 
also allow us to discuss, in section IVIl the implications 
of tracing over spins and study only the entanglement on 
the occupation number degree of freedom compared with 
6]- 

The density matrix in Minkowski coordinates for the 
state (ISTI) is 



P 



M 



p" |0, 0) (0, 0| + pa* |0, 0) (t, tl + |0, 0) (t, i| 



-M7*|o,o)(;,t|- 

-ap* it,t)(t,;i 

-i/3nt,;)(t,;i- 

-l7ni,t)(;,ti- 

-n.d.H.c. 



-/i(5*|o,o)(;,;|- 
-o^i* lt,t)(;,tl 

-h* it,;)(;,ti- 
7-5* i;,t)(;4i- 



"lt,t)(t,tl 

-ar it,t)(;,;i 
-/35*it,;)(;,;i 

(52) 



where n.d.H.c means non-diagonal Hermitian conjugate, 
and represents the Hermitian conjugate only for the non- 
diagonal elements. 

Computing the density matrix, taking into account 
that Rob is constrained to region I of Rindler space-time, 
requires to rewrite Rob's mode in terms of Rindler modes 
and to trace over the unobservable Rindler's region IV. 
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In appendix |^ we compute each term of ((52)) in 
Rindlcr's coordinates and trace over the unobserved re- 
gion IV. Using (jA2p . (jA4p . (jA6| we can easily compute 
the density matrix for Ahce and Rob from (15^ since 
PAR = Tr/y pm, resulting in the long expression 



PAR 



cos" r |0, 0) (0, 0| + sin^ r cos^ r (|0, f) (0, t| 



io,;>(o,;i) + sinV|o,n>(o,ni 



/i COS r 



|o, 0) (t, tl + r|o,o)(t,;i + 7*10,0) a, tl 



-r|o,o)a,;i 



«*|o,;>(t,ni 



II sm r cos r 

io, t) (t, ni +7* io, ;) (;, ni - <5* io, t> a, ni 

+ cos^ r [lap It, t> (t, tl + It, t> (t, tl + "7* 

X It, t> a, tl + It, t> a, ;i + it, ;> (t, ;i 
w It, i) a, tl + P5* It, ;> a, ;i + i7i' It, t> a, ti 

+7ri;,t)a,;i + |5n;,;>a,;i] +sm2r 
■(lap + l/3p) it,t;>(t,t;i + (i7p + l'5p) 
X li, t;> a, t;i + (^7* + it, t;> (;, t;i " 

+n.d.H.c. (53) 



Here the notation is the same than in the r.h.s. of \K2\ : 
|a,r)(a',r'| = \aA) \rR) j {a'j^\ i{r'f(\. Notice that the state, 
which in Minkowski coordinates is pure, gets mixed when 
the observer Rob is accelerated. 

Equation (|53p will be our starting point, from which 
we will study different entanglement settings and how 
Unruh decoherence affects them. 

To begin with we will compute how acceleration affects 
the entanglement of spin Bell states when Alice and Rob 
share a maximally entangled spin state and Rob accel- 
erates. In Minkowski coordinates that means choosing 
specific coefficients in ([5T|) . particularly, for Bell states 
we should choose 



V 2 

|V^±)^/3 = ±7=-i= 



(54) 
(55) 



and the rest of the other coefficients equal to zero. For 
such states in Minkowski coordinates, the density ma- 
trix of Alice and Rob considering that Rob undergoes an 
acceleration a is obtained from (155 



0± 
Par 



cos^ r 



(it,t)(t,ti±it,t)(;,;i±it,t)(t,ti 
i;,;)(;,;i) +sin2r(it,t;)(t,tii 
i;,t;)(;,t;01 (se) 



Par 



cosV( It, ;) (t, ii ± It, t) (t, tl ± It, t) (t, tl 
i;,t)(;,ti) +sin2r(it,tt)(t,t;i 
lt,t;)(;,tt01 (57) 



Notice that, in this case, Alice would have a qubit and 
Rob would have a qutrit, since for his mode he could 
have three different possible orthogonal states: particle 
spin-up, particle spin-down and particle pair. 

To characterize its entanglement we will use the nega- 
tivity [2^ 1 normalized to one (we can multiply it by a 
constant in order to have negativity equal to one for 
a maximally entangled state). Therefore, to have neg- 
ativity equal to 1 for a Bell state we define it as twice 
the addition of all the negative eigenvalues of the partial 
transpose density matrix -which consists on transposing 
Rob's qutrits- 



Par 



cos2r(it,t)(t,ti±it,;)a,ti±i;,t)(;,ti 
it,;)(;,;i) +sin2r(it,t;)(t,tti 
i;,t;)(;,t;01 (58) 



Par. 



cos 



r(it,t)(t,;i±it,t)a,;i±i;,;)(t,ti 

|;,t)(;,tl) +sin2r(|t,t;)(t,ttl 

i;,t;)(;,t;01 (59) 



We can write p'ar'^ matricially in the basis 

{it,t),it,i),ii,t),ii,i),it,ti),i;,t;)} 



/ cos r 




V 





± cos^ r 

± cos^ r 

cos^ 













• 2 

sm r 














• 2 

sm r 



\ 



J 

(60) 

which have the same expression than p^ji'^ in the basis 
{It, i) , It, t) , li, i) , li, t) , It, U) , \i, ti)}. ■ Therefore the 
four Bell states will have the same eigenvalues which are: 



Ai 



A2 — A3 : 
A4 — A5 

As- 



■ cos^ r 



sin^ r 



cos^ r 



(61) 



Since r = arctan (e^'^'S") a — > =^ r — > and a — ?> 
oo=>r— >-7r/4so that Ag is negative for all values of 
the acceleration. This implies, using Peres criterion [25| , 
that the spin Bell states will be always entangled even in 
the limit of infinite acceleration. 
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Figure 2: Negativity and mutual information as a function 
of tiie acceleration of Rob when R and A share a maximally 
entangled state in Minkowski coordinates. Red dashed line 
is mutual information for all the spin Bell states. Black 
dotted line is mutual information for the Minkowski state 
^ (|00) + |t4-)) B.nd blue solid line is negativity for both, Bell 
spin states and ^ (|00) + |t4.>) 



We can readily evaluate the entanglement at the limits 
a — and a — cx) if we compute the negativity (normal- 
ized to one for maximally entangled states), that is to 
say 



\i<0 

Applied to our states we obtain that 
A/'(r) = cos^ r 



(62) 



(63) 



In the limit a — >■ we obtain J\f = 1 which is an expected 
result since a — > is the inertial limit. 

However, in the limit a — )■ oo we obtain A/" = ^, which 
implies that spin entanglement degrades due to the Un- 
ruh effect. Fig. 2 shows the negativity as a function of 
the acceleration of Rob. 

The mutual information, which takes into account 
quantum and classical correlations, is 



Iar = Sa + Sr — Sar 



(64) 



where Sa,r. are the Von Neumann entropies of the partial 
state of Alice and Rob and Sar. is the entropy of the 
whole state. 

For (|56p and (|57|) . the partial states of Alice and Rob 
{pA = Ttrpar, PR = TtaPar) can be expressed matri- 
cially as 



1 



the basis {|t) , 



PR 



1 
1 



3S r 
2 sin^ 



(65) 



(66) 



/ cos- r 



in the basis , ||) , |ti)} for all the Bell states, and 

\ 





PAR 



^ r ± cos^ r 
± cos^ r cos^ r 





sin^ r 



■ 2 

sm r 



the basis {it,t),i;,;),it,n),i;,n)} 



for (j)^ 

and the same expression for ip^ in the basis 
{|t,i),li,t),lt,U),li,ti)}- The entropies of these 
states are 

Sa = 1 

Sr = — cos^ r log2 cos^ r ) — sin^ r log, (sin^ r ) 

2 , 



(68) 



Sar — cos^ r logj (cos^ r) — sin^ r log2 ( - sin^ r 



and the mutual information is 



Iar = 2 cos r 



(69) 



Again we see that in the limit a — mutual information 
goes to 2 and in the limit of infinite acceleration it goes to 
1. The behavior of the mutual information as a function 
of a is shown in fig. 1. 

In it is discussed that Pauli exclusion principle pro- 
tects the on occupation number entanglement from de- 
coherence, and some degree of entanglement is preserved 
even at the limit a — > oo. Here we have obtained a similar 
result for the spin Bell states, showing that spin entan- 
glement is also degraded by Unruh effect. 

Next, we will study the case in which Alice and Rob 
share a different class of maximally entangled state. We 
consider that in Minkowski coordinates we have 



\^) = ±{\0a)\0r) 



It a) Hr)) 



(70) 



which is a maximally entangled state that includes occu- 
pation number entanglement along with spin. We study 
this kind of states as a first analog to the state considered 
in previous literature ([50)1 . This state corresponds to the 
choice 



^ ^ ' = T2 
a = 7 = (5 = 



(71) 
(72) 



in equation (1531) . The density matrix of such a state is 



cos r |0, 0) (0, 0| -I- sin^ r cos^ r (|0, t) (0, tl 

-F|0,4.)(0,;|) +sinV|0,n)(0,ni +cos='r- 
X (|0,0)(t,;| + |t,;)(0,0|) -sin^r cosr 

X (|o, t) (t, ni + It, tt) (0, tl) + cos2 r It, ;) (t, ;i 

+ sin2r|t,ti)(t,t;il (73) 
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Notice the significant difference from tlie Bell spin states; 
considering that Rob accelerates means that, this time, 
Alice has a qubit and Rob has a qu4it. Hence, negativity 
acts only as a measure of distillable entanglement, and 
does not account for the possible bound entanglement the 
system would have [26| . Since in Rindler coordinates the 
state ([75)1 is qualitatively different from the Minkowski 
Bell states ([55]) . (|57l) . it is therefore worthwhile to study 
its entanglement and the mutual information degradation 
as Rob accelerates. 

The partial transpose of ([75)1 cr ~ p^^ is 



cos^ r |0, 0) (0, 0| + sin^ r cos^ r (|0, t) (0, t| 

+ |0,;)(0,;|) +sinV|0,t;)(0,ni +cos3r 
X (|0,i)(t,0| + |t,0)(0,;|) -sin^r cosr 

X (|o, n> (t, ti + It, t> (0, ni) + cos^ r It, ;> (t, ;i 

+ sin2r|t,n>(t,nil (74) 



which is an 8 X 8 matrix, a is diagonal by blocks with 
eigenvalues 



in the basis {|0) , |t) , |i) , Iti)}- The eigenvalues of the 
whole system 6x6 matrix par are 



Ai 
A3 



A2 -0 
1 

2 
1 



■ 2 2 
sm r cos r 



■ 4 
— sm r 
2 



A4 

A5 = -cos^r(l + cos 
Ae 



i sin^ r (1 + cos^ r) 



(79) 



In this case the mutual information as a function of a is 
not proportional to the negativity. Hence, it is different 
from the Bell states (ISBl) . ([57]) . As it can be seen in Fig. 
1 the value of mutual information for (l56l) , (|57)) and (|73)) 
coincide at the limits a — > 0, a — > 00, but are different in 
between, obtaining that /l^"^''" > 7^^'^'=^'=". 

To conclude this section we stress that the same results 
wiU be obtained if the state |t,4.) in ([70)) is replaced by 
any other 1 particle bipartite spin state |s, s'). 



Ai 
A2 
A3 
A4 
As.e 

A7,8 



cos^ r 



2 • 2 

cos r sm r 



. 2 

sm r 



cos^ r 



(75) 



- ^sin^ r cos^ r ± \/ sin"* r cos* r + 4 cos^ 



sin"* r ± V sin^ r + 4 sin"' r cos^ r 



As we can see, As is non-positive and Ag is negative for 
all values of a, therefore the state will always preserve 
some degree of distillable entanglement. If we calculate 
the negativity we will obtain 



7V(r) = cos^ r 



(76) 



which means that for this case, distillable entanglement 
behaves equally than in the previous case, and negativity 
on fig. 1 is equally valid for this state. 

Finally we compute the mutual information of the sys- 
tem whose partial matrices are expressed as 



1 n 
1 



(77) 



in the basis {|0) , |t)}. pn. 



( cos"* r 

sin^ r cos^ r 

V 



cos^ r (sin^ r 



sm r sm r 



+ 1)7 

(78) 



VI. OCCUPATION NUMBER 
ENTANGLEMENT WITH AN ACCELERATED 
PARTNER AND SPIN 1/2 FERMIONS 

The previous work Q on occupation number entan- 
glement between accelerated partners ignored the spin 
structure of the Dirac field modes. It is not possible to 
straightforwardly translate a state like (|5T|) into mere oc- 
cupation number states. This comes about because for 
a state like (|5ip the bipartite vacuum component does 
not have individual spin degrees of freedom as the other 
components do. In other words, by including the vacuum 
state in the superposition (|51l) . the Hilbert space ceases 
to be factorable in terms of individual spin times particle 
occupation number subspaces. 

On the other hand, the bipartite vacuum is a well de- 
fined total spin singlet. Hence, the Hilbert space is fac- 
torable with respect to the total spin of the system A — R 
and the occupation number subspaces. Accordingly, to 
reduce the spin information in the general density ma- 
trix (15^ we will be forced to consider a factorization of 
the Hilbert space as the product of the total spin and 
occupation number subspaces. 

If we do such a factorization we could consider that we 
are not able to access to the information of the total spin 
of the system A—R and then, we should trace over total 
spin degree of freedom. 

The equivalence between the standard basis (occupa- 
tion number-individual spin) and the new basis (occupa- 
tion number-total spin) is given^ in equations (j80p and 



^ The pair state in the same mode can only be a singlet of total 
spin due to anticommutation relations of fermionic fields 
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(.pIr) 



|0,0) = |00)|5) 


io,n> = 


= |02> |^> 


|0,t) = |01)|i?+) 


|o,;) = 


101) \D_) 


|t,0) = |10)|i?+) 


i;,o) = 


|10)|^-) 


|t,n) = |12)|^?+ 


) i;,n) = 


\12)\D^ 


lt,t> = |ll>|r+> 


i;4> = 


|11>|T_> 




11) [l^o) + 


\S)] 




11) [\To) - 


\S)] 



(80) 



(81) 



where we are using the basis |na «{,) | J, J^) and the 
triplets, doublets and the singlet are denoted as 



\T+) = 


\J = 


l.Jz = 


1) 


\T-) = 


\J = 


l,Jz = 


-1> 


\To) = 


\J = 


l,Jz = 


0) 


D+) = 


\J = 


1/2, J. 


-1/2) 


D-) = 


\J = 


1/2, J. 


= -1/2) 


\S) = 


\J = 


0, J. = 


0) 



(82) 

If we rewrite the general state ([5T|) in this basis we obtain 
I*) 



M|00)|5)+a|ll)|T+) + ^|ll)|ro) 



V2 



P-1 
V2 



ll)\S)+5\n)\T^) 



(83) 



And the general state when Rob is accelerated (|53t in 
terms of this new basis after reducing the information on 
the total spin by tracing over this degree of freedom is 



Par = {J,Jz\par\J,Jz) 



(84) 



Which results in a state in the occupation number basis 
whose entanglement decoherence could be studied and 
compared with the results in reference Q in which spin 
is ignored: 



Par 



cos^ r |00) (00| + 2 sin^ r cos^ r |01) (01| 
sin" r 1 02) (02 1 



+ ^cos^r ( ^* S^* |00)(11| 



V2 



/3-7 



V2 

-sin2r|12)(12| 



|11)(00| + (1 -/r) cos^r|ll)(ll| 



(85) 



we can readily compute the partial transpose ct" 



^i' cosV|00)(00| + 2sin2r cos2r|01)(01| 
sin"* r 1 02) (02 1 



^cos^r ( ^ ^ |01)(10| 



n/2 



/3-7 



V2 

+ sin2r|12)(12| 
whose eigenvalues are 

\ 2 4 

Ai = /I cos r 



|10)(01| + (1 -M^) cos^r |11)(11| 



(86) 



2 • 4 

/I sm r 



A3 = (1 - /i^) cos^ r 
A4 = (l-^^)sin^r 



^^5,6 = cos^ r I /i^ sin^ r ± /ia/ /i^ sin" r + cos^ r^— 



all the eigenvalues are non- negative except Ae < 0. The 
negativity ((62)) is, in this case. 



Af = 2 cos^ 



/i^ sin^ r — /iy //^ gij^^ _|_ qqs^ r 



1/3 -7l^ 



(87) 

which depends on the proportion of singlet \ j3 — 7|/\/2 of 
the 1 11) component in the state (|5T|). When there is no 
singlet component (/3 = 7) the negativity is zero. Indeed 
in the limit a — > (Minkowskian limit) 



AAo = y2|/i||/?-7| 



(88) 



That shows that the maximally entangled Minkowski oc- 
cupation number state (Negativity = 1) arises tracing 
over total spin when the starting state is 



|vi/) = i=|o,o)±i[it,;)-i;,t)] 

or, in the occupation number-total spin bases 
1 



I*) 



V2 



m\s)±\ii)\s) 



(89) 



(90) 



That means that, for occupation number entanglement, 
the only way to have an entangled state of the bipartite 
vacuum |00) and the one particle state |11) of a Dirac 
field is through the singlet component of total spin for 
the |11) component. 

On the contrary, the state 



I*) 



1 r 
71 



|00)|5)±|ll)|To,±) 



(91) 



will become separable after tracing over total spin due to 
the orthonormality of the basis ([50]) . ([5T|) . 

We have established that the Minkowski maximally en- 
tangled state for occupation number arises after tracing 
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Figure 3: Negativity (blue solid line) and mutual information 
(red dashed line) as a function of the acceleration of Rob when 
R and A share an occupation number maximally entangled 
state (|89|l in Minkowski coordinates after tracing over total 
spin 



over total spin in a state as Now we will compute 

the limit of the negativity when the acceleration goes to 
c» in order to see its Unruh decoherence and to compare 
it with the results for occupation number entanglement 
from 0. 

Taking a->oo=»r^7r/4in dST]) 



1 



m'-Vm^ + m' 1/3-71' 



(92) 



Therefore, for the maximally Minkowski entangled state 
we have ^ = l/-\/2, |/3 — 7I = 1 and the negativity in the 
limit is 



J^oo = 



(93) 



This result shows that when we are reducing the total 
spin information, looking at the occupation number en- 
tanglement alone, we see that it is more degraded by 
Unruh effect than when we considered spin Bell states in 
the previous section. More importantly, the occupation 
number entanglement is more degraded than in Q , where 
the spin structure of the modes was considered nonexis- 
tent. This happens because considering spin structure 
of each mode, occupation number 2 is allowed. Hence, 
Pauli exclusion principle protection of the entanglement 
is weaker than in [§1 where the spin is not considered. 
The negativity dependence on the acceleration is shown 
in fig. 3 

We can also compute the mutual information for the 
state (|55t as we did in the rest of the cases. Its analyt- 
ical expression is quite long and has no special interest, 
but we can see the dependence of Iar for the Minkowski 
maximally entangled state (|89l) with the acceleration in 
fig. 3, obtaining that /^^ = 2 and = 1/2 



VII. CONCLUSIONS AND COMMENTS 

It is known 0, Q that Unruh decoherence degrades en- 
tanglement of occupation number states of fields. Here 
we have shown a richer casuistic that appears when we 
take into account that each Dirac mode has spin struc- 
ture. This fact enables us to study interesting effects 
(such as Unruh decoherence for spin Bell states) and de- 
velop new procedures to erase spin information from the 
system in order to study occupation number entangle- 
ment. 

Along this work we have analyzed how a maximally 
entangled spin Bell state losses entanglement when one 
of the partners accelerates. We have seen that, while in 
Minkowski coordinates Alice and Rob have qubits, when 
Rob accelerates the system becomes a non-pure state of 
a qubit for Alice and a qutrit for Rob. In this case spin 
entanglement for a Dirac field is degraded when Rob ac- 
celerates. However some degree of entanglement survives 
even at the limit a 00. 

A first analog to the well studied state (l/-\/2)(|00) -I- 
|11)) but including spin could be, for instance, 
(l/^/2)(|00) + |tl>). This state, unlike the deceivingly 
similar spin Bell states, becomes a qubit xqu4it when 
Rob accelerates. Nevertheless, distillable entanglement 
degrades in the same way as for spin Bell states. 

We have also introduced a procedure to consistently 
erase spin information from our setting preserving the 
occupation number information. We have done it by 
tracing over total spin. The maximally entangled oc- 
cupation number state is obtained from the total spin 
singlet after tracing over total spin. Finally we have 
shown that its entanglement and mutual information is 
more degraded than in Q where the spin structure of 
Dirac modes was neglected. A reasonable physical ar- 
gument for this result is that, in our setting, occupation 
number 2 is allowed for the Dirac field modes, and hence, 
there is a broader margin for entanglement degradation 
by Unruh effect. 

The thermal noise (|48l) is obtained when dealing with 
a two dimensional space-time and massless fields. Mass 
gap and transverse degrees of freedom modify the count- 
ing statistics that is no longer given by thermal noise, but 
replaced by the so-called Rindler noise 22|, that the the 
space-time dimension. In this work we were concerned 
with the specific issues associated to the spin degree of 
freedom, so the restriction to massless fields in two di- 
mensional space-time adopted here, as well as the single 
mode approximation, allows a direct comparison with the 
previous works j^, 161] which considered massless spinless 
fields in 2D under those approximations. 

As a matter of fact, having more than 2 space-temporal 
dimensions and having massive fields may introduce rel- 
evant effects. Allowing the possibility of having mo- 
mentum off the acceleration direction and having mas- 
sive fields we would obtain a spread of Minkowski modes 
over Rindler frequencies p3 |. If we carry out the single- 
mode approximation, the spread of Minkowski modes 
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into Rindler modes can be neglected even for higher di- 
mensions [l|, 0, [l^ , but if we want to relax such an ap- 
proximation (as for the discussion in the next paragraph), 
those effects should be considered in order to account 
for the entanglement in non-inertial frames. It would be 
worthwhile to study those effects in future articles. 

Another very interesting point that deserves further 
study is the fact that when we consider more than one 
populated mode of the complete Minkowski vacuum ^ 
instead of the single mode approximation, the margin for 
Unruh degradation increases as we could have in princi- 
ple a larger number of levels that can be excited by the 
Thermal/Rindler noise. One could think that these cases 
would be quite similar to the bosonic case Q where the 
margin for Unruh decoherence is so broad that no en- 
tanglement survives at the limit of a — > oo. Something 
similar would apply as well if we relax the single mode 
approximation allowing a small spread in both Rob's de- 
tector and populated levels such that we consider a con- 
tinuum of accessible levels. These topics will be inspira- 
tion for future works. 
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Appendix A: Some Minkowski operators expressed 
in Rindler's region I 

To help with the calculations of the density matrix 
associated with ([51]) it is useful to compute firstly the 
trace over IV on all the operators that compose ([5^. 
Using equation ([^^ we have 

|0,0)(0,0| = |0a) [cos^r |0)jO)j^ -f e'-^sinr cosr 
x(|t>,i;)/v + i;>/lt)/v) + e^'*sinV 

xin)/in>/y]®H.c. (Ai) 



tracing over IV: 

Tr/y |0,0)(0,0| = cos^r |0,0)(0,0| 

+ sin2rcos2r(|0,t>(0,t| + |0,i)(0,i|) 
+ sin4r |0,n)(0,n| (A2) 

where notation is different in each side of the equality: 
bras and kets in l.h.s. are referred to Minkowski coor- 
dinates for Alice and Rob \s,s') = \sa) \s'j^) and in the 
r.h.s. they are referred to Alice's mode in Minkowski co- 
ordinates and Rob's mode in Rindler's region I \s, s') = 

\sa) \s'ii)j. 

In the same way, using expressions (UU , ([33]) we have 

|0,0)(s,s'| = \0A)[cos^r\0)j\0)jy + e"''smrco8r 
X(|t>/l;>/y + l;>/lt)/y) + e2^'^sin2r 
x|n)/ln>/y] {s\{cosri{s'\iv{0\ 
+£'e-''*' sin ri{n\iv{s'\) (A3) 

where e' = 1 if s =t and e' = — 1 if s Now, tracing 
over IV 

Tr/v |0,0)(s,s'| = cosV |0,0)(s,s'| -Ksin^r cos7■ 
x(<5,,t|0,;)(.s,n|-<5.'JO,t)(.s,nl) 

(A4) 

notation here is the same than in (jA2l) . 
Again, using expression PH) we get 

|si,S2)(s3,S4| = \si)[cosr \s2)j\0)jy+S2e~"''smr 
X \U)i\s2)iv] (ssl [cosr/(s4|/v(0| 
+e4e-*'^sinr,(n|7y(s4|] (A5) 

and tracing over IV gives 

Tr/y |si, S2)(S3, S4I = COS^ r |si, S2)(S3, S4I 

+6,^s^ sin2r|si,n)(s3,til 

(A6) 

again, notation here is the same than in (jA2[) . 
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